Spatial optical solitons supported by mutual focusing 
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We study composite spatial optical solitons supported by two- wave mutual focusing 
induced by cross-phase modulation in Kerr-like nonlinear media. We find the families 
of both single- and two-hump solitons and discuss their properties and stability. We 
also reveal remarkable similarities between recently predicted holographic solitons in 
photorefractive media and parametric solitons in quadratic nonlinear crystals. 
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One of the standard physical mechanisms support- 
ing the existence of spatial optical solitons is the beam 
self-focusing due to nonlinearity-induced change of the 
medium refractive index (see, e.g., Refi and references 
therein). In such a case, an optical beam induces an ef- 
fective waveguide and becomes self-guided in a nonlinear 
optical medium. 

A rather different physical mechanism is responsible 
for the so-called mutual beam focusing which can be ob- 
served, for example, in photorefractive two- wave mixing 2 . 
In this case, a single beam diffracts because its self- 
induced focusing is weak or negligible. However, the same 
beam can demonstrate focusing in the presence of the 
other beam, due to the mutual cross-phase-modulation 
interaction in a Kerr-like medium with the third-order 
nonlinear susceptibility. As was predicted recently, the 
effect of mutual focusing can be responsible for the forma- 
tion of a new kind of optical spatial soliton, the so-called 
holographic solitor^. Such a soliton is created in the ab- 
sence of the self-action effect, solely due to the beam 
cross-coupling through the induced Bragg reflection. Be- 
cause the holographic soliton is described by two modes 
in the frames of the coupled-mode theory, it can be re- 
garded as a special case of vector solitons, well studied 
in nonlinear optics. 

The concept of a vector soliton was originally associ- 
ated with nonlinear interaction of two polarizations in a 
birefringent Kerr medium. In this case, both self- and 
cross-phase modulation of the beams produce a bound 
state of two orthogonally polarized field components in 
the form of a self-trapped two-component vector solitor*. 
The existence of spatial vector solitons can be understood 
through the concept of induced waveguiding. Indeed, 
when one of the soliton components creates an effective 
waveguide in a nonlinear medium via the self-modulation 
effect, the second component can be localized as a guided 
mode of this waveguide, which creates a composite state 
in the form of a vector soliton even beyond the bifurcation 
point 5 . This concept allows to introduce, in a simple and 
straightforward way, both single- and multi-hump vec- 
tor solitons, once it is accepted that the soliton-induced 
waveguide can support higher-order modesi. 

In the case of vector solitons due to mutual focusing 
the physics of beam coupling cannot be grasped by the 



simple waveguiding argument. In particular, it is not 
intuitively clear what would be the properties of such 
solitons when the powers of the two constituents are 
significantly dissimilar. Since the holographic solitons 
have been found only in the special case when two opti- 
cal beams have equal amplitudes (and the corresponding 
mathematical problem reduces to a scalar equation), we 
wonder if the mutual self-focusing may generate stable 
composite solitons in a Kerr-like medium in a broad range 
of experimental parameters, and if they can be compared 
with other types of solitons studied earlier. 

First of all, we note that mutual focusing is also known 
to support parametric spatial solitons in quadratic non- 
linear medieQ. In general, such three-wave paramet- 
ric solitons consist of two coupled low-frequency beams 
and the sum-frequency component. In this case, the 
phase-matching energy-dependent interaction between 
the beams is a key physical mechanism for supporting 
the parametric solitons. One of the waves can be gener- 
ated inside the crystal, however at least two waves at the 
input are required to generate a three-wave parametric 
soliton. This feature suggests that there may be a strong 
similarity between different types of solitons supported 
by mutual focusing. 

In this Letter we analyze the composite spatial solitons 
supported by mutual beam focusing in a Kerr-like nonlin- 
ear medium in the absence of the self-action effects. We 
predict the existence of continuous families of single- and 
two-hump composite solitons, and discuss their stability 
and interaction. 

In order to describe the mutual focusing of two optical 
beams in a nonlinear medium with the third-order sus- 
ceptibility, we consider the model for holographic solitons 
recently introduced by Cohen et al£, 



.du 
dz 
.dv 

oz 



a u 
dx 2 
d 2 v 
dx 2 



\v\ 2 u 



2 _|_ |„|2 

u\ 2 v 



1 + \uy + \v 

|2 



1 



= 0, 



= 0, 



(1) 



where u(x, z) and v(x, z) are the normalized envelopes of 
the two interacting waves in a waveguide geometry. The 
waves are coupled only through the effective nonlinear 
cross-phase modulation terms and the saturation effect 



is also taken into account^. Since the mode coupling is 
effectively incoherent in this approximation, the partial 
mode powers, P\ = f_°° \u\ 2 dx and P 2 = /J 1 " 00 |i>| 2 dx, 
are conserved independently. 

At low powers, nonlinear interaction is weak and the 
diffraction-induced spreading of both beams is observed. 
However, at higher powers mutual focusing of the beams 
becomes important. Spatial solitons correspond to a 
balance between nonlinear and diffraction effects; the 
corresponding solutions of Eq. are sought in the 
form: u(x,z) = U(x) exp(iaiz + icpi) and v(x,z) — 
V(x) exjp(ia 2 z + if > 2 ), where cti t2 are the nonlinear prop- 
agation constants and 1^1,2 are the arbitrary phases. The 
real amplitude profiles U(x), V(x) satisfy the following 
system of coupled equations, 
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We look for localized solutions describing bright soli- 
tons for which light is localized in the central region, 
(U, V)(x — > ±00) — ► 0. In a special degenerate case 
when ati = a%, the profiles of two components coincide, 
U(x) = V(x), and such solutions have been recently de- 
scribed by Cohen et al£. In this Letter, we study a gen- 
eral case when the amplitude profiles of two waves do 
not coincide. Generally speaking, the two-wave soliton 
solutions can be found only numerically, but the exis- 
tence region of symmetric single-hump solitons can be 
evaluated analytically. Indeed, from the symmetry prop- 
erty for single-hump solitons, {U, V}(x) = {U, V}(— x), 
it follows that the point x — is a maximum, and we 
use Eq. J5J) to obtain the corresponding inequalities for 
the soliton peak amplitudes U m — U(0) and V m = V(0), 
V 2 (l + U 2 n + V 2 )-\ > ai and Ul (1 + U 2 m + > a 2 . 

These two inequalities are compatible only when ai. 2 > 
and (ai + a 2 ) < 1, so that these conditions define the ex- 
istence region of the single-hump vector solitons. The 
soliton peak amplitudes are bounded from below, 

m 1 - (at + a 2 ) ' m 1 - (ai + a 2 ) ' W 

Therefore, the peak intensity of one wave is proportional 
to the propagation constant of the other component; this 
relation reflects the fact that the two- wave soliton is sup- 
ported solely by the mutual focusing effect. Both inten- 
sities rapidly increase and the saturation effect becomes 
important close to the boundary of the existence region 
where (a\ + a 2 ) — > 1. 

To find the soliton families, we solve Eq. J5J numer- 
ically and find localized solutions everywhere inside the 
analytically estimated existence region on the parame- 
ter plane (ai,a 2 ) shown in Fig. ^top) as a triangular 
domain in the inset. The characteristic dependencies of 
the soliton partial powers are presented in Fig.^top) as 
functions of the propagation constant a 2 . The character- 
istic examples of the two- wave soliton profiles are shown 




FIG. 1: Top: Powers Pi and P2 of the soliton components as 
functions of the soliton parameter ai at fixed ot\ — 0.1. Inset 
shows the existence domain on the parameter plane, the solid 
line ai = ai corresponds to the previously discussed solitons^. 
Thin horizontal line indicates the bistable solutions. Bottom: 
Examples A, B, and C of the two-component single- hump soli- 
tons supported by mutual focusing. Solid and dashed curves 
correspond to the fields v and u, respectively. 
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FIG. 2: (a)-(d) Examples of two-component two-hump soli- 
tons supported by mutual focusing for different vales of a.2 
(indicated on the plots) and fixed ai = 0.5. 



in Figs, ^bottom), and they correspond to the points 
A, B, and C marked on the power plot. In agreement 
with the analytical results J3J), the amplitude of the u 
component increases for larger values of the propagation 
constant a 2 , and the corresponding power dependence 
P\{a 2 ) is monotonic. On the other hand, the power de- 
pendence ^2(012) for the v component exhibits bistabil- 
ity. The power increases for a 2 — > because in the low- 
saturation regime the amplitude of the v component is 
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FIG. 3: Collision of two in-phase vector solitons in (a) high 
(ai = 0.39, « 2 = 0.59) and (b) low (ai = 0.039, a 2 = 0.059) 
saturation regimes. 

almost fixed by the value of a%, whereas the beam width 
grows as y/o2. In spite of such a non-monotonic depen- 
dence of the power component, the two-wave solitons cor- 
responding to both branches are stable. The reason is the 
following. In the system under consideration two partial 
powers are conserved independently, as the mutual inter- 
action is incoherent. In this case, the soliton stability is 
defined by the generalized Vakhitov-Kolokolov criterion 
for multi-component solitons^ which, as we have verified 
for our model, predicts soliton stability. Thus, single- 
hump vector solitons supported by mutual focusing are 
stable. 

We now search for two-hump vector solitons for which 
the components u and v have symmetric unipolar and an- 
tisymmetric bipolar profiles, respectively. We find that 
such solitons appear only when a,\ < cti- This is exactly 
one half of the existence region for one-hump vector soli- 
tons in the inset shown in Fig.^top). The characteristic 
examples of two-hump solitons are presented in Fig. |3 
where we show the modification of the soliton profiles 
when one of the parameters is fixed. As the boundary of 
the existence domain is approached, i.e. «2 - * ot\, the 
separation between two humps increases, and the two- 
hump solutions can be considered as a bound state of 
two single-hump vector solitons. However, for small ai 
the unipolar component becomes much wider, and it car- 
ries more power than the bipolar one. In this limit the 
structure of two-hump vector solitons resembles that of 
parametric solitons in nonlinear quadratic media™. 

Our results suggest that there should exist a link be- 
tween holographic and three- wave parametric solitons. 
Let us consider the properties of parametric solitons 6 - 9 
when the phase mismatch Afc is large. Then, the gener- 
ated sum-frequency wave has a small amplitude, which 
can be found using the cascading approximation, Ej, ~ 
EiE 2 /Ak. After substituting this expression into the 



governing equations for the low-frequency beams, we ob- 
tain the model which is mathematically equivalent to 
Eqs. QJ in the low-saturation limit. Therefore, many 
similarities between holographic and parametric solitons 
are not a coincidence, but rather a manifestation of the 
similar physical origin of such solitons - the effect of mu- 
tual focusing. 

We would like to stress that, despite some similari- 
ties, the soliton dynamics in the low- and high-saturation 
regimes can be substantially different, as is clearly illus- 
trated in Figs. 01a,,b) for the soliton collisions. In both 
cases solitons merge, however only in the high-saturation 
regime strong oscillations develop, see Fig. EI a). This 
happens because, in the regions with high beam inten- 
sity, the medium response is almost linear, and persis- 
tent beating between several modes can develop. On the 
other hand, in the low-saturation case oscillations quickly 
decay due to strong radiation, see Fig. Etb). 

In conclusion, we have analyzed two-component spa- 
tial solitons supported by mutual focusing of two optical 
beams in the absence of the self-phase modulation ef- 
fects. One of the possible realization of these solitons 
is through the holographic focusing effect recently dis- 
cussed by Cohen et a/£. We have demonstrated the exis- 
tence of the continuous families of single- and two-hump 
vector solitons, and have discussed their stability. Ad- 
ditionally, we have emphasized a key difference between 
the formation of the convectional vector solitons and the 
two-component solitons supported by mutual focusing, 
but also revealed their links to the parametric quadratic 
solitons. 
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